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ABSTRACT 

A  combined  experimental  and  numerical  study  has  been  carried  out  to  examine  the 
influence  of  thickness  on  the  surface  and  through-thickness  stress  distributions  of 
aluminium  alloy  specimens  with  central  holes.  Thermographic  and  strain-gauge 
measurements  were  made  on  several  thicknesses  of  specimens,  the  central  holes  of 
which  were  either  reamed  or  cold-expanded  by  3.6%.  The  experimental  results  were 
compared  with  those  from  numerical  analyses  of  the  specimens  by  three-dimensional 
finite-element  modelling. 
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Effect  of  Plate  Thickness  on  the  In-Plane  and 
Through-Thickness  Stresses  at  a  Hole 


Executive  Summary 


Holes  exist  in  metallic  aircraft  structures  for  either  design  or  construction  reasons  and 
their  stress  concentrating  effect  is  such  that  they  are  common  regions  of  fatigue  crack 
initiation. 

Theoretical  studies  have  indicated  that  the  stress  concentration  factor  (maximum  local 
stress/ remote  stress)  of  a  hole  varies  through  the  thickness  of  a  plate  and  is  sensitive  to 
the  value  of  hole  diameter  to  plate  thickness  ratio  (D/t),  and  to  Poisson's  ratio.  Theory 
has  also  shown  that  the  assumption  of  plane  stress  is  adequate  for  thin  plates  which 
have  a  thickness  less  than  one  quarter  the  hole  diameter.  For  thick  plates -which  have 
a  thickness  of  at  least  twice  the  diameter  of  the  hole -the  three-dimensional  stress 
concentration  factor  should  be  used. 

To  assess  the  validity  of  the  various  solutions,  an  experimental  study,  including 
thermoelastic  and  strain-gauge  measurements,  has  been  carried  out  at  AMRL.  This 
report  presents  the  results  of  the  study  which  examined  the  influence  of  plate 
thickness  on  the  three-dimensional  stress  distribution  of  aluminium  alloy  specimens 
with  central  holes.  The  variations  in  surface  stresses  and  through-tihe-thickness 
stresses,  with  specimens  of  thickness  10,  20  or  35  mm,  have  been  investigated.  Three- 
dimensional  elastic  finite-element  analyses  were  also  conducted. 

The  experimental  data  compared  well  with  the  numerical  finite-element  work  for  the 
surface  of  the  specimens.  The  corresponding  through-thickness  data  were  not  as 
clearly  in  accordance,  due  to  experimental  problems.  The  numerical  finite-element 
work  undertaken  in  this  report  validates  the  trends  of  other  published  finite-element 
and  analytical  studies,  that  as  the  thickness  of  a  plate  with  a  hole  increases,  the 
maximum  stress  concentration  factor,  which  is  in  the  middle  of  the  plate,  increases  and 
the  minimum  value  of  stress  concentration  factor,  which  is  at  the  surface  of  the  plate, 
decreases.  It  is  shown  that  caution  must  be  taken  in  engineering  calculations  for  thick 
plates  with  holes  — the  common  use  of  a  two-dimensional  handbook  stress 
concentration  factor  may  be  too  conservative  at  the  surface,  but  unconservative  at 
mid-thickness. 
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Notation 

D 

Hole  diameter  (mm) 

E 

Yotmg's  modulus  (MPa) 

H 

Specimen  height  (mm) 

Kt 

Elastic  stress  concentration  factor  =  (oy/ Oo) 

W 

Specimen  width  (mm) 

k 

Calibration  factor 

r 

Hole  radius  (mm) 

t 

Specimen  thickness  (mm) 

£x/  Sy/  Sz 

Normal  strain  in  the  x,  y  and  z  directions,  respectively 

V 

Poisson's  ratio 

Obulk 

Bulk  stress  (ox  +  Oy  +  Oz)  (MPa) 

Oo 

AppHed  gross-area  stress  (MPa) 

CJx/  CTy,  Oz 

Normal  stress  in  the  x,  y  and  z  directions,  respectively  (MPa) 

Abbreviations 


CE  Cold-Expanded 

FAST  Focal-plane  Array  for  Synchronous  Thermography 

NCE  Non-Cold-Expanded 

PAFEC  Program  for  Automatic  Finite  Element  Calculations 

SCF  Stress  Concentration  Factor 

SD  Standard  Deviation 

SPATE  Stress  Pattern  Analysis  by  the  measurement  of  Thermal  Emission 
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1.  Introduction 

Recent  work  relevant  to  the  assessment  of  structural  integrity  of  aircraft  has 
considered  the  influence  of  three-dimensional  stresses  at  holes.  Several  papers  have 
presented  analytical  solutions  for  the  three-dimensional  stress  field  at  a  hole  in  a  large 
plate  of  arbitrary  thickness;  [1-5]  being  a  selection.  Historical  summaries  of  this  subject 
can  be  found  in  papers  by  Folias  and  Wang  [1],  and  Youngdahl  and  Sternberg  [2].  In 
the  analytical  papers  the  plate  is  assumed  to  be  of  a  homogeneous,  isotropic,  elastic 
material,  to  be  infinite  in  width  and  height,  and  have  a  remote  tensile  stress  applied. 
These  studies  have  shown  that  the  stress  concentration  factor  varies  through  the 
thickness  of  the  plate  and  is  sensitive  to  the  ratio  of  hole  diameter  to  plate  thickness 
(D/t)  and  Poisson's  ratio  (v),  and  that  the  Yoimg's  modulus  of  the  material  has  no 
effect  on  the  value  of  the  factor.  Other  findings  are;  that  the  assumption  of  plane  stress 
is  adequate  for  ratios  of  D/t>4,  and  for  ratios  of  D/t<0.5  the  three-dimensional  stress 
concentration  factor  should  be  used. 

A  finite-element  study  on  this  problem  has  been  presented  by  Shivakumar  and 
Newman  [6].  As  in  the  analytical  studies,  the  plate  was  assumed  to  be  of  a 
homogeneous,  isotropic,  elastic  material,  however,  ihe  plate  was  finite  in  width  and 
height.  Several  load  cases  were  covered  (tension,  bending  and  pin-loading),  and  both 
straight  and  countersunk  holes  were  considered.  (It  should  be  noted  that  in  the  course 
of  this  present  study,  an  error  was  discovered  in  this  paper.  The  figures  and  text 
mention  that  the  plate  (with  a  straight-shank  hole)  has  a  height  and  width  of  5  times 
the  diameter  of  the  hole,  whereas  the  height  and  width  are  equal  to  10  times  the  hole 
diameter.  The  authors  have  confirmed  this  error  and  a  subsequent  ESDU  data  sheet 
(93030)  based  on  [6]  is  to  be  amended.) 

These  several  studies  [1-6]  do  not  appear  to  be  associated  with  any  confirmatory 
experimental  program.  To  validate  the  various  numerical  solutions,  an  experimental 
study  has  now  been  carried  out  at  AMRL  which  examines  the  influence  of  plate 
thickness  on  the  three-dimensional  stress  distribution  of  aluminium  alloy  specimens 
with  central  holes.  The  variations  in  surface  stresses  and  through-the-thickness 
stresses,  with  specimens  of  thickness  10,  20  or  35  mm,  have  been  investigated.  The 
stress  distribution  in  the  specimens,  which  each  had  a  reamed-only  or  a  cold-expanded 
central  hole  of  diameter  22.6  mm,  was  obtained  via  thermographic  measurements 
taken  by  a  system  called  FAST  (Focal-plane  Array  for  Synchronous  Thermography) 
[7].  Some  of  the  specimens  also  had  strain-gauge  measurements  taken.  Three- 
dimensional  elastic  finite-element  analyses  were  conducted  for  calibration  of  the 
thermographic  data,  and  for  comparison  with  the  strain-gauge  data  and  previous 
analytical  work. 
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2.  Theoretical  Analyses 

The  main  consideration  of  the  analytical  papers  is  the  variation  in  stress  concentration 
factor  (oy/oo  =  Kt)  through  the  thickness  of  the  specimen.  However  the  variation  in 
(Oz/Oo)  through  the  thickness  of  the  specimen  is  also  required  to  enable  comparison 
with,  and  calibration  of,  the  experimental  FAST  data  which  are  output  in  digital  levels 
that  range  from  -128  to  127.  Three  of  the  papers  referenced  above  have  determined  the 
Oz  component,  however  they  did  not  consider  two  of  the  important  conditions  in  the 
present  examination,  namely,  a  plate  and  hole  of  finite  dimensions,  and  v  =  0.33  (for 
aluminium).  The  strain-gauge  locations  for  the  experimental  stress  determinations 
were  unsuitable  for  calibrating  FAST,  hence  finite-element  analyses  were  undertaken 
with  the  added  benefit  of  determining  stresses  for  the  exact  material  and  geometry 
tmder  consideration.  These  three-dimensional  elastic  finite-element  (FE)  stress 
analyses  were  undertaken  using  the  PAFEC  suite  of  programs,  level  8.1,  run  on  a 
HP9000-755  computer. 

2.1  Comparison  with  Published  Results  for  a  Circular  Hole  in  a  Large 
Plate 

Two  additional  FE  analyses  were  undertaken  to  give  confidence  in  the  modelling 
method  used;  one  to  compare  with  the  corrected  numerical  study  of  [6],  and  the  other 
to  compare  with  Folias  and  Wang's  analytical  work  on  an  infinite  plate  [1]. 

2.1.1  Numerical  Study 

The  geometry  and  material  properties  of  the  PAFEC  model  were  identical  to  one  of  the 
cases  studied  in  [6].  The  plate  which  was  under  a  tensile  load,  had  a  height  and  width 
of  lOD,  that  is,  ten  times  the  hole  diameter  (refer  to  Fig.  1  for  axis  and  dimension 
notation)  and  v  =  0.3,  and  a  similar  mesh  to  that  in  [6]  was  used.  The  PAFEC  stress 
concentration  factors  (Kt)  for  the  selected  case  of  D/t=2  were  between  0.28%  and  0.65% 
lower  than  those  reported  in  [6].  Several  verifications  of  the  FE  analysis  were 
completed  and  it  was  concluded  that  the  model  was  suitably  accurate.  As  an 
independent  check,  the  same  case  was  repeated  with  the  NASTRAN  FE  package; 
version  68.1.0,  November  15  1994,  run  on  an  IBM  RISC  System-6000,  [8].  A  similar 
mesh  to  the  PAFEC  model  was  used  and  the  Kt  values  for  the  comparison  with  [6] 
were  less  than  0.15%  higher  than  the  corresponding  PAFEC  results,  refer  to  Table  1. 


Table  1.  Comparison  of  values  ofSCF  as  determined  by  various  FE  programs  for  D/t=2  and 
W=H=10D 


z/t 

Stress  Concentration  Factor 

1  %  Difference  | 

Ref.  [6] 

PAFEC 

NASTRAN 

Ref.  [6]  & 
PAFEC 

Ref.  [6]  & 
NASTRAN 

PAFEC  & 
NASTRAN 

0 

3.2 

3.191 

3.195 

0.281 

0.156 

0.125 

0.5 

2.94 

2.921 

2.925 

0.646 

0.510 

0.137 
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(a)  View  of  specimen  surface 
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(b)  Section  through  y-z  plane 


Figure  1 .  Axis  and  notation  for  circular  hole  in  a  plate 

2.1.2  Analytical  Study 

The  infinite  plate  case  examined  analytically  by  Folias  and  Wang  where  D/t  =  2  was 
modelled  by  the  same  large  plate  geometry  of  section  2.1.1  (ie,  height  and  width  equal 
to  lOD)  witii  V  =  0.33.  The  PAFEC  Kt  values  for  this  case  (of  D/t=  2)  were  on  average 
2.7%  higher  than  those  reported  in  [1].  Again  an  independent  check  was  made  by 
repeating  the  comparison  using  the  NASTRAN  FE  package.  NASTRAN  gave  Kt  values 
which  were  less  than  0.14%  higher  than  the  corresponding  PAFEC  results  (refer  to 
Table  2)  hence,  any  difference  between  the  PAFEC  and  NASTRAN  models  is 
negligible.  Part  of  the  discrepancy  between  the  analytical  and  numerical  "infinite" 
plate  could  be  attributed  to  the  finite  width  effect  (assessed  as  ^1.3%  from  [9]).  The 
balance  of  the  difference  («»1.4%)  can  readily  be  attributed  to  the  finite  height, 
numerical  error  of  the  FE  package,  the  accuracy  of  [1],  or  a  combination  of  the  three. 


Table  2.  Comparison  of  values  ofSCF  as  determined  by  analytical  and  FE  models  for  D/t=2, 
and  V^-H-IOD  for  FE  models 


z/t  1 

Stress  Concentration  Factor 

j  %  Difference  | 

Ref.  [1] 

PAFEC 

NASTRAN 

i  Ref.  [1]& 

1  PAFEC 

Ref.  [1]  & 
NASTRAN 

PAFEC  & 
NASTRAN 

0 

3.1 

3.191 

3.195 

2.935 

3.065 

0.125 

0.5 

!  2.85 

2.921 

2.925 

i  2.491 

2.632 

0.137 

3 


DSTO-TR-0330 


2.2  Numerical  Analysis  of  Test  Specimens 

The  three  non-cold-expanded  experimental  specimens.  Fig.  2,  were  modelled  using  20- 
noded  isoparametric  solid-brick  elements.  The  mesh  of  the  model  of  the  10  mm  thick 
specimen  is  shown  in  Fig.  3  and  is  comprised  of  2513  elements.  The  PAFEC  data  file 
for  this  model  is  listed  in  Appendix  A.  Due  to  symmetry,  only  one-eighth  of  the 
structure  was  modelled.  The  meshes  for  the  20  mm  and  the  35  mm  thick  specimens 
were  comprised  of  4308  and  4185  elements,  respectively.  Appropriate  boundary 
condition  constraints  were  placed  at  all  planes  of  symmetry.  A  remote  stress  of  30  MPa 
was  applied  to  the  x-z  plane  at  a  location  corresponding  to  the  placement  of  the 
fatigue  machine  grips  on  the  specimen,  as  schematically  depicted  in  Fig.  2.  The  data 
obtained  from  these  analyses  included  bulk-stress  values  on  the  surface  and  down  the 
bore  of  the  hole  as  appropriate  to  compare  with  experimental  readings. 


Region  of 
specimeii 
in  grips 


specimen  thickness 
t=10,  20  or  35 


Figure  2.  Geometry  of  specimens  (Dimensions  in  mm) 
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3.  Experimental  Analyses 


3.1  Specimen  Details 

The  six  specimens  used  in  this  investigation  were  manufactured  froin  45  thick 

rolled  plates  of  2214-T651  aluminium  alloy  and  consisted  of  a  rectangular 

centrally  located  hole.  The  specimen  configuration  is  shown  m  2/  and  the 

speSme^ns  have  constant  thicknesses  of  either  10,  20  or  35  mm  (two  of  each  thickness). 

The  central  hole  in  half  of  the  specimens,  one  of  each  thickness,  was 

fCE)  by  3  6%  (using  the  Boeing  split-sleeve  system  as  discussed  m  [10])  fte 

lightly  mamed  to  22.6  mm  diameter.  The  remainder  had  a  reamed-only  hole  (NC  ) 

22.6  mm  in  diameter. 

To  experimentally  determine  the  through-thickness  and  surface  stresses  two  metho^ 
wereLployed,Lmely,  strain-gauge  measurements,  and  bulk-stress 
via  a  thainographic  system.  Both  sets  of  measuremei^s  were  made  whUe  the 
specimen  was  loaded  in  a  standard  servohydraulic  test  machme. 

3.2  Strain-Gauge  Stress  Analysis 

3.2.1  Testing  Details 

The  four  10  mm  and  35  mm  thick  specimens  were  strain-gauged  down  the  bo-e  of Jhe 
hole  to  measure  y  and  z  strains  with  120  ohm  strip-gauges  havmg  a  pitch  of  2.03 
The  35  mm  specimens  were  also  strain-gauged  on  the  surface  to  measure  x  and  y 
strains  as  per  the  schematic  of  Fig.  4.  The  x  and  z  direction  gauges  had  an  active  ^d 
wth'of  0  7874  mm  and  width  of  0.8128  mm,  and  the  y  direction  gauges  had  an  achve 
grid  length  of  0,7874  mm  and  width  of  1.778  mm.  Strain  readings  were  taken  every 
flvlPa  ^ling  from  zero  to  30  MPa  to  -30  MPa  to  zero  (these  applied  shesses  were 
based  on  gross  area).  Not  all  of  the  gauges  were  functioning  correctly,  therefore  the 
data  gained  were  limited.  The  data  acquisition  system  used  was  a  digital  volt-meter 

type  3455A  HP. 

3.2.2  Calibration 

As  mentioned  previously,  a  schematic  of  the  location  of  the  gauges  is  shown  in  Fig  4. 
Due  to  the  small  difference  in  the  positioning  of  the  x  and  y  strip  gauges  on  the 
surface,  and  also  the  y  and  z  strip  gauges  down  the  bore  of  the  hole,  mterpolation  o 
the  strain  readings  was  performed.  Using  these  interpolated  values,  the  stresses  were 
calculated  using  the  following  equations; 


on  the  surface; 


ax= 


E(Sx  +  V6y) 


ay= 


E(8y  +  V8x) 

1-v' 
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down  the  bore  of  the  hole; 


ay  = 


E(£y  +  V£z) 

l-v" 


(5z  = 


E(£z+  VSy) 

1-v' 


where  E  =  73,084  MPa  and  v  =  0.33 

It  should  be  noted  that  a  strain-gauge  reading  is  an  average  over  the  area  of  the  grid  of 
the  gauge  though  it  is  usually  taken,  as  at  present,  to  be  the  value  at  the  centre  of  the 
grid.  Hence  the  stresses  calculated  as  above,  are  also  average  values. 


10  mm  thick  specimens 
Hole 

2  strain-gauges  in  y  direction 

3  strain-gauges  in  z  direction 


35  mm  thick  specimen 
Surface 

10  strain-gauges  in  x  direction 
10  strain-gauges  in  y  direction 

Hole 

7  strain-gauges  in  y  direction 
7  strain-gauges  in  z  direction 


FIGURE  4.  Schematic  of  strain-gauge  locations 
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3.3  Thermoelastic  Stress  Analysis 
3.3.1  Testing  details 

Thermographic  measurements  were  taken  of  the  bore  of  the  hole  and  of  the  face  of  all 
six  specimens  (on  the  face  opposite  to  the  strain-gauges  where  appropriate).  This  was 
performed  by  the  FAST  system  with  the  fatigue  test  machine  applying  an  alternating 
stress  of  ±30  MPa  at  a  frequency  of  5  Hz.  The  resultant  scans  were  486  by  510  pixels  in 
size.  The  system  used  in  this  experimental  analysis  is  based  on  a  real-time  digital  video 
processor  and  a  focal-plane  array  thermal  imager  [7],  which  detects  the  minute 
changes  in  infra-red  radiation  caused  by  the  thermoelastic  effect.  Changes  in 
temperature  of  less  than  0.01°  C  can  be  measured,  with  the  specimens  painted  matt 
black  to  improve  their  emissivity.  This  system  is  able  to  deliver  high  resolution  fuU- 
field  thermoelastic  stress  scans  of  the  specimen  with  a  reduction  in  observation  time  of 
between  one  and  two  orders  of  magnitude  [11]  compared  to  the  single  sensor 
thermographic  systems  (such  as  SPATE)  which  create  comparable  stress  plots  by 
scanning  one  point  at  a  time. 


3.3.2  Calibration 

As  mentioned  previously,  the  thermographic  measurements  for  each  specimen 
consisted  of  a  surface  scan  and  a  scan  down  the  bore  of  the  hole.  Two  horizontal  lines 
of  data  were  extracted  from  the  surface  scans;  one  corresponding  to  the  location  of  the 
X  and  y  surface  gauges  (but  on  the  opposite  face),  and  the  other,  at  location  y  =  45  mm. 
A  line  of  data  through  the  thickness  of  the  specimen  parallel  to  the  z  axis, 
corresponding  to  the  location  of  the  y  and  z  hole  strain-gauges  (but  on  the  opposite 
side)  was  obtained  along  the  bore  of  the  hole  from  the  other  scan.  These  line  scans 
were  adjusted  for  aberration. 

3.3.2.1  Surface  Data 

Various  boundary  conditions  to  model  the  mechanics  of  the  fatigue-machine  grips 
were  implemented  into  the  FE  runs.  From  these  results  it  was  ascertained  that  within  a 
vertical  distance  of  45  mm  from  the  centre  of  the  hole,  the  bulk  stress  distribution  in 
the  plate  was  unaffected  by  the  grips.  On  the  surface  of  the  specimens  the  bulk  stress 
was  equal  to  (ox  +  Oy),  as  Oz  =  0.  The  plots  of  bulk  stress  for  the  numerical  models 
(Fig.  5 -which  also  includes  the  two-dimensional  result  for  comparison)  indicate  that 
at  this  distance  (of  45  mm)  thickness  effects  are  small.  Thus  the  calibration  location  was 
selected  as  line  y  =  45  mm.  The  spatial-average  numerical  bulk-stress  value  at 
y  =  45  mm  was  calculated  for  the  specimen  models  of  different  thicknesses.  The 
average  uncalibrated  bulk-stress  value  from  the  FAST  data  was  then  determined  from 
the  horizontal  line  scans  at  the  calibration  location  for  each  specimen.  The  calibration 
factor  (k)  for  each  specimen  was  thus  calculated  as  the  ratio  of  the  averages  of 
numerical  FE  bulk  stress  to  uncalibrated  experimental  bulk  stress  (refer  to  Table  3). 
The  FAST  data  for  each  specimen  were  then  multiplied  by  the  corresponding  factor  k. 
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Figure  5.  Bulk-stresses  parallel  to  the  y  axis  at  x=0,  determined  by  FE,for  specimens  of 
different  thickness 


Table  3.  Determination  of  calibration  factors  for  FAST  surface  bulk  stresses 


Specimen 

FAST  uncalibrated 
average  stress 

FE  spatial-average 
stress 

Factor 

(k) 

10  mm,  NCE 

47.56 

30.21 

0.635 

10  mm,  CE 

41.88 

30.21* 

0.721 

20  mm,  NCE 

45.49 

30.20 

0.664 

20  mm,  CE 

41.23 

30.20* 

0.733 

35  mm,  NCE 

50.23 

30.24 

0.602 

35  mm,  CE 

44.50 

30.24* 

0.680 

*  These  values  are  from  a  finite-element  model  of  a  NCE  specimen 
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3.3.2.2  Through-Thickness  Data 

The  imcalibrated  bulk-stress  data  obtained  from  the  through-thickness  scans  for  each 
specimen  were  noisier  than  the  surface  data  and  a  quadratic  curve-fit  was  applied  to 
each  Une  scan.  At  this  location  down  the  bore  of  the  hole  the  bulk  stress  is  equal  to 
(oy  +  Oz)  as  Ox  =  0.  The  caHbration  factor  (k)  was  calculated  by  determining  the  spatial- 
average  of  the  bulk-stress  along  a  line  in  the  z  direction  (at  y  =  0)  as  obtained  by  FE 
analysis,  and  dividing  by  the  average  of  the  corresponding  curve-fitted  uncahbrated 
FAST  data  (Table  4).  The  FAST  data  for  each  specimen  were  then  multiplied  by  the 
corresponding  factor  k. 


Table  4.  Determination  of  calibration  factors  for  FAST  through-thickness  bulk  stresses 


Specimen 

FAST  imcalibrated 
average  stress 

FE  spatial-average 
stress 

Factor 

(k) 

10  mm,  NCE 

101.0 

99.46 

0.985 

10  mm,  CE 

109.3 

99.46* 

0.910 

20  mm,  NCE 

70.47 

104.4 

1.481 

20  mm,  CE 

56.40 

104.4* 

1.851 

35  mm,  NCE 

76.19 

109.6 

1.439 

35  mm,  CE 

76.55 

109.6* 

1.432 

*  These  values  are  from  a  finite-element  model  of  a  NCE  specimen 


4.  Results  and  Discussion 

The  majority  of  the  experimental  and  numerical  FE  results  are  summarised  in  Figs.  6  to 
11.  Figures  6,  7  and  8  are  normalized  plots  of  the  numerical  and  experimental  surface 
data  for  the  10,  20  and  35  mm  thick  specimens  respectively,  where  the  bulk  stresses  are 
normalized  by  the  applied  stress.  Figures  9,  10  and  11  are  normalized  plots  of  the 
numerical  and  experimental  through-thickness  data  for  the  10,  20  and  35  mm  thick 
specimens  respectively;  the  FAST  data  are  represented  by  the  associated  curve-fit  and 
standard  deviation  (SD).  As  with  the  surface  data,  the  bulk  stresses  in  Figs.  9-11  are 
normalized  by  the  applied  stress. 

As  mentioned  in  Section  2.2,  the  numerical  models  are  of  the  non-cold-expanded 
specimens  only,  therefore  the  FE  curve  in  the  cold-expanded  plots  (Figs,  (b)  of 
Figs.  6-11)  is  of  the  non-cold-expanded  (NCE)  model  of  corresponding  thickness.  It 
should  be  noted  that  in  the  absence  of  a  full  elastic/plastic  three-dimensional  FE 
analysis  of  a  cold-expanded  hole  and  relevant  strain-gauge  data,  calibration  of  the 
FAST  data  is  complicated.  Therefore,  the  results  presented  for  the  cold-expanded 
specimens  show  the  trends  of  the  data,  not  absolute  values.  These  specimens  were 
tested  to  experimentally  investigate  the  effect  of  cold-expansion  on  the  surface  and 
through-thickness  stresses  at  a  hole. 
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4.1  Comparison  of  Experimental  and  Finite-Element  Surface  Data 

From  the  surface  lines  of  data  in  the  x  direction  at  y  =  0  (Figs.  6-8)  it  can  be  seen  that 
the  experimental  FAST  data  agree  well  with  the  numerical  FE  curves  (in  most  cases). 
At  the  edges  of  the  hole  and  plate,  the  FAST  data  points  drop  rapidly  away  from  the 
numerical  curves  due  to  the  nature  of  the  FAST  system;  the  data  points  are  averaged 
by  those  around  them.  Hence,  a  sharp  thermal  representation  at  an  edge  cannot  be 
attained.  For  the  10  and  20  mm  thick  specimens,  the  FAST  data  indicate  that  the  bulk 
stress  is  higher  in  the  NCE  specimens.  The  FAST  data  of  the  35  mm  thick  specimens 
tend  to  indicate  the  opposite,  however,  the  strain-gauge  readings  closest  to  the  hole  for 
these  specimens  (Fig.  8)  indicate  that  the  bulk  stress  is  lower  for  the  cold-expanded 
specimen,  as  for  the  10  and  20  mm  thick  specimens.  However,  due  to  a  lack  in 
consistency  between  the  FE  and  FAST  NCE  data,  there  is  not  much  confidence  in 
drawing  a  conclusion  from  the  FAST  CE  data.  Apart  from  a  few  points,  the  strain- 
gauge  measurements  (available  for  the  35  mm  thick  specimens  only)  tend  to  follow  the 
data  of  the  FE  model  and  the  FAST  data. 


4.2  Comparison  of  Experimental  and  Finite-Element  Through- 
Thickness  Data 

The  experimental  data  gained  from  the  through-thickness  bulk-stress  measurements 
are  of  poorer  quality  than  desired.  As  mentioned  earUer,  the  FAST  data  were  noisy  and 
had  to  be  represented  by  a  curve  of  best-fit.  The  SD  which  is  indicated  in  Figs.  9-11  is 
quite  large  for  the  35  mm  thick  specimens.  On  the  whole,  the  general  trends  between 
the  numerical  and  experimental  data  for  the  different  thickness  specimens  are  in 
agreement,  ie,  the  maximum  value  of  bulk  stress  tends  to  occur  in  the  middle  of  the 
specimen  (z/t  =  0)  and  the  minimum  on  the  surface  (z/t  =  +0.5),  however,  the 
correlation  between  the  magnitudes  of  the  numerical  and  experimental  data  is  a  Httle 
unsatisfactory.  Due  to  the  nature  of  calibration  it  is  not  possible  to  make  a  profoimd 
statement  regarding  the  comparison  of  the  NCE  and  CE  results. 

4.3  General  Trends  of  Finite-Element  Through-Thickness  Results 

It  can  be  seen  from  the  numerical  FE  curves  of  Figs.  6-8  that  as  the  thickness  of  the 
specimen  increases,  the  maximum  value  of  normalized  bulk  stress,  which  occurs  at  the 
edge  of  the  hole,  reduces.  This  is  hence  mirrored  in  the  values  of  the  normalized  bulk 
stress  at  the  surface  in  the  corresponding  through-thickness  curves  of  Figs.  9-11.  These 
figures  also  show  that  the  maximum  value  of  the  normalized  bulk  stress  occurs  in  the 
middle  of  the  plate,  increasing  with  increase  in  thickness.  The  Oz  component  of  this 
bulk  stress  accounts  for  a  small  percentage  of  the  total,  the  maximum  value  ranging 
from  13%  to  56%  of  the  remote  applied  stress  (oo),  whereas  the  maximum  value  of  the 
CTy  component  ranges  from  331%  to  338%  of  Oo.  For  the  geometry  of  specimens  tested 
the  maximum  Oz  and  Cy  occurred  in  the  middle  of  the  plate.  The  maximum  Oz  value 
increased  with  increase  in  thickness,  however,  the  maximum  Oy  value  occurred  in  the 
20  mm  thick  specimen. 
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The  through-thickness  FE  data  generated  in  the  present  work  are  summarised  in 
Fig.  12,  and  are  expressed  as  stress  concentration  factor  versus  normalized  distance 
through  the  thickness.  Two  additional  curves  are  plotted  on  this  figure;  an  extra  FE 
run  with  D/t=  8,  and  theoretical  two-dimensional  data  from  Peterson  [9],  both  for  the  r 
and  W  dimensions  of  the  specimens  tested.  As  mentioned  in  Section  2.1,  the  numerical 
data  of  Fig.  12  are  slightly  higher  than  that  of  the  published  work,  however,  the  trends 
are  similar.  From  Fig.  12  it  can  be  seen  that  for  the  large  values  of  D/t  the  gradient  of 
the  curves  is  small  near  the  surface  of  the  specimen,  and  the  maximum  values  of  Kt 
occur  in  the  centre  of  the  specimen.  The  middle  value  of  D/t  (1.13)  has  a  relatively 
higher  value  of  maximum  Kt  and  the  gradient  of  the  curve  near  the  surface  is 
increased.  The  curve  for  the  thickest  specimen  (ie,  lowest  D/t)  has  the  largest  gradient 
near  the  surface  of  the  specimen,  however,  not  the  highest  maximum  Kt.  Therefore,  as 
the  value  of  D/t  is  increased,  the  maximum  value  of  Kt  in  the  centre  of  the  specimen 
and  the  gradient  of  the  curve  near  the  surface  of  the  specimen  both  reduce,  tending 
towards  a  straight  line,  ie,  the  plane  stress  case.  The  variation  in  von  Mises  stress  down 
the  bore  of  the  hole  with  ratio  of  D/t  shows  the  same  trend  as  the  Kt  curves,  however,  it 
is  more  pronounced.  Thus,  as  the  thickness  of  the  specimens  increase,  the  maximum 
value  of  von  Mises  stress  moves  away  from  the  centre  of  the  plate  towards  the  surface, 
which  perhaps  is  an  indicator  of  the  location  of  fracture,  as  mentioned  in  [1].  Findings 
from  theoretical  papers  also  indicate  that  for  ratios  of  D/t>  4  the  assumption  of  plane 
stress  is  adequate,  and  the  transition  region  from  plane  stress  to  plane  strain  occurs 
between  D/t  ratios  of  0.5  to  4.  Though  not  covering  a  full  range  of  D/t  ratios,  the 
present  numerical  FE  work  confirms  these  conclusions. 


Normalized  bulk  stress  Normalized  bulk  stress 
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(a)  Non-cold-expanded  specimen 


(b)  Cold-expanded  specimen 

Figure  7.  Comparison  of  experimental  and  numerical  normalized  bulk-stresses  parallel  to  the 
X  axis  at  y=0  for  20  mm  thick  specimens 


DSTO-TR-0330 


.  FE  results 
Strain-gauge  data 
.  Fit  to  FAST  data 
SD  =  0.56 


Normalized  distance  along  bore  of  hole  (z/t) 


(a)  Non-cold-expanded  specimen 


.  FE  results 
Strain-gauge  data 
.  Fit  to  FAST  data 
SD  =  0.58 


-0,3  -0.2  -0.1  0  0.1  0.2 
Normalized  distance  along  bore  of  hole  (z/t) 


(b)  Cold-expanded  specimen 


Figure  11.  Comparison  of  experimental  and  numerical  normalized  bulk-stresses  along  the  bore 
of  tlm  hole  for  35  mm  thick  specimens 
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5.  Conclusions 


1.  The  elastic  three-dimensional  finite-element  work  undertaken  validates  the  trends 
of  the  published  FE  and  analytical  studies,  that  as  the  thickness  of  a  plate  with  a 
hole  increases,  the  maximum  stress  concentration  factor,  which  is  in  the  middle  of 
the  plate,  increases  and  the  minimum  value  of  stress  concentration  factor,  which  is 
at  the  surface  of  the  plate,  decreases.  Thus,  care  must  be  taken  in  engineering 
calculations  for  thick  plates  with  holes— the  use  of  a  two-dimensional  handbook 
stress  concentration  factor  may  be  too  conservative  at  the  surface,  but 
unconservative  at  mid-thickness. 

2  The  three-dimensional  finite-element  analyses  have  indicated  that  Oz  accounts  for 
only  a  small  percentage  of  the  bulk  stress  through  the  thickness  of  the  specimen 
compared  to  Oy.  The  maximum  value  of  Oz  ranged  from  13%  to  56  %  of  the  remote 
applied  stress,  whereas  the  maximum  value  of  the  Oy  component  ranged  from 
331%  to  338%  of  the  remote  applied  stress.  For  all  thicknesses  of  specimen  tested 
(10,  20  and  35  mm)  the  maximum  of  Oz  and  Oy  occurred  in  the  middle  of  the 
plate— the  maximum  Oz  value  increased  with  increase  in  thickness,  however,  the 
maximum  Oy  value  occurred  in  the  20  mm  thick  specimen. 

3.  The  experimental  FAST  data  compared  well  with  the  numerical  FE  work  for  the 
surface  of  the  specimens.  The  corresponding  through-thickness  data  were  not  as 
clearly  in  accordance.  Apart  from  a  few  points,  the  strain-gauge  trends  appeared 
to  agree  with  the  FE  models  and  the  thermoelastic  data.  However,  it  would  seem 
that  for  better  experimental  results  for  the  through-thickness  work,  either  smaller 
and  more  reUable  strain-gauges  and/or  FAST  measurements  with  reduced  noise 

are  needed. 
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Appendix  A 

PAFEC  data  file  for  10  mm  thick  specimen 

C  PAFEC  3D  GEOMETRY  FOR  A  HOLE  IN  A  PLATE 
C  MODEL  IS  1/8  SYMMETRIC,  MATERIAL:  AL  ALLOY  2214-T651 
C  THICKNESS:  10  mm,  HOLE  DIAMETER:  22.6  mm 
C 

CONTROL 

FULL . CONTROL 

PHASE=1 

PHASE=2 

PHASE=4 

PHASE=6 

PHASE=7 

PHASE=9 

CONTROL . END 

C 

NODES 


NODE. NUMBER 

X 

Y 

Z 

1 

0 

11.3 

0 

2 

4.324322786 

10.43983872 

0 

3 

7 . 990306628 

7 . 990306628 

0 

4 

10 . 4398387 

4.324322786 

0 

5 

11 . 3 

0 

0 

6 

20 . 0 

0 

0 

7 

48.0 

0 

0 

8 

0 

20.0 

0 

9 

20.0 

20 . 0 

0 

10 

48.0 

20.0 

0 

11 

0 

115.0 

0 

12 

48.0 

115.0 

0 

13 

20.0 

115.0 

0 

14 

0 

0 

0 

15 

0 

11.3 

5.0 

16 

4.324322786 

10.43983872 

5 . 0 

17 

7 . 990306628 

7 . 990306628 

5.0 

18 

10.4398387 

4.324322786 

5.0 

19 

11.3 

0 

5.0 

20 

20.0 

0 

5 . 0 

21 

48.0 

0 

5 . 0 

22 

0 

20.0 

5.0 

23 

20.0 

20.0 

5.0 

24 

48.0 

20 . 0 

5 . 0 

25 

0 

115.0 

5 . 0 

26 

48 . 0 

115.0 

5.0 

27 

20.0 

115 . 0 

5 . 0 

28 

0 

0 

5 . 0 

C  END  OF  NODES 
C 

C  BRICK  CONFIGURATION 

PAFBLOCKS 

PROPERTIES=ll 

ELEM=37110 

TYPE=1 


BLOCK 

GROUP 

N1 

N2 

N5 

TOPOLOGY 

1 

1 

1 

4 

6 

1,3,15,17,8,9,22,23,2,0,0,16 

2 

2 

2 

4 

6 

3,5,17,19,9,6,23,20,4,0,0,18 

3 

3 

3 

7 

6 

6,7,20,21,9,10,23,24 

4 

4 

1 

5 

6 

8,9,22,23,11,13,25,27 

5 

C 

C  END 

5  3 

OF  PAFBLOCKS 

5 

6 

9,10,23,24,13,12,27,26 
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C 

MESH 

REFERENCE 

1 


SPACING. LIST 
5 

1,1, 1,1, 1,1, 1,1,1, .5, .5 

1,1, 2, 2, 4, 4 

.45, .45,1,1,2,2,2,2 

2,3,5,5,5,5,5,5,5,5,10,10,10,10,10 
0.4, 0.4, 0.6, 0.6, 1,1,1 
.5, .5, 1,1, 1,1, 1,1, 1,1,1 


C  END  OF  MESH 
C 

MATERIAL 
MATERIAL . NUMBER 
11 
C 

C  END  OF  MATERIAL 
C 

RESTRAINTS 

NODE . NUMBER  PLANE 

1  1 

5  2 

15  3 

C 

C  END  OF  RESTRAINTS 
C 

SURFACE . FOR . PRESSURE 
PRESSURE . VALUE  NODE 
-30.0  11 

C 

C  END  OF  PRESSURE 


E 

73084.5 


DIRECTION 

1 

2 

3 


PLANE 

2 


C 

END . OF . DATA 
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